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An exact stochastic model for the thermalisation of quantum states is 
proposed. The model has various physically appealing properties. The dy- 
qq . namics are characterised by an underlying Schrodinger evolution, together 

0^ | with a nonlinear term driving the system towards an asymptotic equilibrium 

state and a stochastic term reflecting fluctuations. There are two free param- 
eters, one of which can be identified with the heat bath temperature, while 
the other determines the characteristic time scale for thermalisation. Exact 
expressions are derived for the evolutionary dynamics of the system energy, 
the system entropy, and the associated density operator. 
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I. INTRODUCTION 



The construction of a microscopic description for dynamical systems out of equilibrium 
is a topic of considerable interest in statistical mechanics. While many studies have been 
pursued in this direction M, the subject remains refractory owing to its nonlinear character. 



In particular, a realistic framework must be simple enough to be tractable, yet sufficiently 
sophisticated to capture some of the physical essentials. Our goal here is to propose in this 
q-( spirit an elementary but compelling stochastic model that characterises the thermalisation 

process for a quantum system in an arbitrary given initial state. 

The approach here, though carried out largely within the framework of standard quan- 
tum theory, involves a reexamination of the traditional hypotheses forming the basis of 
quantum statistical mechanics. In particular, we shall take the view that a quantum sta- 
tistical ensemble, or 'state', is correctly described not simply by a density matrix p^, but 
rather by a probability distribution on the space of pure quantum mechanical states. For 
measurements associated with ordinary quantum mechanical observables this description 
reduces to the usual one, for which the unconditional expectation of a linear observable is 
given by the familiar trace formula involving p%. The advantage of the present formulation, 
however, is that we can give meaning to the concept of ensembles and other mixed states in 
nonlinear quantum mechanics @J^]. Furthermore, the fact that the space of pure states has 
a symplectic structure allows us to introduce the concept of a quantum phase space T, and 
hence the use of probability distributions on F as ensembles. In particular, for thermal equi- 
librium in certain situations we propose the distribution given in (|T2"D below, which we call 
the canonical T-ensemble. The corresponding density matrix differs in some respects from 
the conventional density matrix of quantum statistical mechanics, and thus the applicability 
of the theory may in principle be testable. 
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Our first step is to derive the Fokker-Planck equation for the time-dependent probability 
density function associated with a general Brownian motion with drift on the space of pure 
quantum states. Since T is compact, we can apply a theorem of Zeeman to show that there 
is a special class of drift terms for which, given any initial condition, the resulting asymp- 
totic stationary solution is the canonical T-ensemble described above. The nonequilibrium 
dynamics we consider here is governed by a nonlinear Schrodinger evolution coupled to a 
thermal noise term, described by the stochastic differential equation (H). This choice is the 
simplest possible that leads to thermal equilibrium for arbitrary initial states. 

In what follows we use the dynamical law (|TJ]) to derive a formula for the total system 
energy (|I8D, leading to a formula for the heat capacity (|19|). These ideas are applied, by 
way of illustration, to the case of a spin one-half particle in a heat bath. Then we derive an 
expression for the dynamics of the system entropy, given in ([24]), from which we are able to 
deduce a fundamental thermodynamic inequality (pop- Finally we study the evolution of the 
density matrix associated with the underlying process. We show that the time development 
of the density matrix involves a higher moment of the projection operator onto pure states. 



II. QUANTUM HAMILTONIAN DYNAMICS 

We begin by considering the quantum mechanical phase space CP n [0. This is the space 
of pure quantum states, given by rays through the origin in the associated Hilbert space. 
The quantum phase space can be viewed as a 2n-dimensional real manifold T, endowed with 
the Fubini-Study metric g a b- The reason that T plays the role of a quantum phase space is 
that Hamilton's equations on T can be lifted to the Hilbert space to give the Schrodinger 
dynamics. This follows from the fact that V has a natural symplectic structure, given by a 
nondegenerate antisymmetric tensor Q a b, compatible with the metric g a b in the sense that 
Va^bc = and Q a b^ bc = S^, where Q ab = g ac g bd Q c d and V a denotes the standard covariant 
derivative, satisfying V a gbc = 0. Then the Schrodinger dynamics take the Hamiltonian form 

dx a = 2Q ab V b Hdt , (1) 

where dx a /dt is tangent to the quantum phase space trajectory. The Hamiltonian function 
H(x) on r is given by the expectation of the energy operator at each point x G T. More 
specifically, the expectation of the Hermitian operator with respect to the Hilbert space 
vector ip a (a,/3 = 0, 1, • • • ,n) is given by H(x) = H^ipaip 13 /V> 7 V> 7 • Since each point x G T 
corresponds to an equivalence class {\ip a , A G C — 0} for some Hilbert space vector ip a , it 
follows that H(x) is defined globally on T. Conversely, one can show that such functions 
correspond to global solutions of the equation 

V 2 H = (n + l)(H-H) , (2) 

where V 2 is the Laplace-Beltrami operator on T. Here H = H^/(n + 1), being the trace 
of the Hamiltonian. Thus H is the uniform average of the energy eigenvalues, corresponding 
to the thermal equilibrium energy of the system at high temperature. 
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III. STOCHASTIC DYNAMICS 



Now we generalise the Schrddinger dynamics by consideration of a diffusion process Xt 
taking values in T. For this purpose we make use of standard techniques of stochastic 
differential geometry [fH|. We shall examine the case of a Brownian motion with drift on 
r, determined by the covariant stochastic differential equation 

dx a = p a dt + KofdWf , (3) 

where k is a constant, fi a is a vector field, and the vectors cr" (i — 1, 2, - • • , 2n) constitute an 
orthonormal basis in the tangent space of T such that g ab = afa b 5 %: > and (T^a b g a b = S^. Here 
dx a is the covariant Ito differential ||, and the standard 2n-dimensional Wiener process W[ 
satisfies dW\dWl = 5^dt. In local coordinates x a (a = 1, 2, • • • , 2n), the Ito differential is 
given by 

dx a = 5 a a {dx a + ^K 2 T a c g hc dt) . (4) 

Here 5% is a coordinate basis for the given patch of T, and S a the dual basis, such that for 
the covariant derivative V a £ fe of a vector field £ a with components £ a = 5 a £, a we can write 

df a 

te(v 6 r) = ^ + rSe^ . (5) 

The Ito differential (|]) is constructed in such a way that the stochastic differential equation 
(0) is fully tensorial. 

Suppose 0(x) is a smooth function on T, and define the process f = 4>{x t ). It follows 
from Ito's lemma that d<p t = V a 0cfe a + \ V a V b<t>dx a dx b . Then the relation dx a dx b = n 2 g ab dt 
implies 

d4> t = (/i a V a + -k 2 V 2 0)cft + nS/a^dWt , (6) 

where V 2 = g afc V a Vb0. Since the process Xt is Markovian, there is a normalised density 
function p(x, t) for the state at time t, characterised by a partial differential equation subject 
to specified initial conditions. The expectation of the process <j) t is thus 



E[<f> t ] = / <j>(x)p(x,t)dV , (7) 
J r 

where dV denotes the volume element on r. On the other hand, it follows from Ito's lemma 
(|) that 

& = 0o + f (V V a + ^ 2 V 2 ^) <is + k ay a <f) dWl . (8) 

Here the integrands are valued at time s at the point x s G T. Since the second integral 
above has vanishing expectation, we obtain 



E\6 t ] = <b Q + E 



j\p a V a <P+^K 2 S/ 2 <P)ds] . (9) 
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Hence differentiating ([?]) and @ with respect to t and equating the results, we have 

2 

J p(x, t)4>{x)dV = | r (/x a V a + y V 2 0)p(x, t)dV, (10) 

where p = dp/dt. Integrating by parts, and using the fact that the resulting relation must 
hold for all 4>(x), we find that p(x,t) satisfies 

^p(x,t) = -W a (p a p) + 1 -K 2 V 2 P . (11) 

This is the covariant Fokker-Planck equation for the density function on F, corresponding 
to the stochastic process ([3]). The solution of flTID thus characterises the distribution of the 
diffusion (|3|) at time t, given an initial distribution p(x,0). In the case of a singular initial 
distribution, e.g., an initially pure state, we interpret p(x,t) as a generalised function, i.e., 
it has to satisfy ( |10"D for all smooth <p(x). 

The foregoing results are valid on arbitrary compact Riemannian manifolds. We are 
concerned, however, with the case where F is the quantum phase space CP n , endowed with 
the Fubini-Study metric. Since F is compact, equation (|TT|) may admit nontrivial asymptotic 
(t — > oo) stationary solutions. We shall examine a simple situation in which this is the case, 
namely, when the drift is given by a gradient flow p a = — \K 2 f3V a H generated by H(x), 
where (3 is a parameter. Then we can use a theorem of Zeeman to show that there is a 
unique stationary solution for ([TT|) of the form 

p(x) = «P(-y» (12) 

where = J r exp(— /3H(x))dV. It follows that the probability density function 

p(E) for the distribution of the energy function H(x) is given by the formula p(E) = 
Z- 1 (l3)n(E)exp(-(3E), where 

Q(E) = J r 6(H(x) - E)dV (13) 

is the density of states in F for which E < H(x) < E + dE, and Z{0) is the partition 
function. Another way of characterising the distribution ([12]) is that it maximises the entropy 
S = — J r p(x) In p(x)dV for a given value of system energy U = J r p(x)H(x)dV. The 
theorem of Zeeman implies in the present context that, (0) is the asymptotic distribution of 
the process Xt under essentially arbitrary initial conditions, and that p(E) is the asymptotic 
energy distribution. 

We have so far considered a process of the form ([3]) with the drift 2 f3V a H. This 

process, as such, does not yet make reference to the Schrodinger dynamics. If we view the 
drift term as a nonlinear correction to the underlying quantum evolution, we can represent 
the complete dynamics according to the prescription 

dx a = (2Q ab V b H - ^K 2 (3V a H^j dt + na^dW l t , (14) 

which in the limit k — > reduces to the ordinary Schrodinger dynamics. Due to the anti- 
symmetry of Q a b, inclusion of the symplectic term does not affect the resulting asymptotic 
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state fljjD, since p(x) is a function of H(x). Hence, the analysis above shows that, given 
an arbitrary initial state (pure or general) on the quantum phase space, the dynamical law 
(ITJJ) necessarily takes that state into the thermal equilibrium state ([12]). In particular, the 
process (|14|) is ergodic on the energy surface, and asymptotically leads to a uniform dis- 
tribution on each such surface. We note that ( |H| ) involves two parameters, (3 and k. The 
stationary solution for p(x, t) depends only on (3, which we identify as the inverse of the 
heat bath temperature. The parameter k, which has the units s _1//2 , and may depend on (3, 
determines the thermalisation time scale. 



IV. HAMILTON! AN PROCESS 

Our next step is to study the stochastic process H t = H(x t ) associated with the Hamil- 
tonian function. From Ito's lemma (|6|) for the process fll4|) we find 

dH t = ^K 2 (V 2 H t - 0V t )dt + K\/ a HafdWi . (15) 

Here V t = V(x t ), where V(x) = g ab V a HVbH is the squared energy uncertainty at the point 



x G r, conditional on the pure state to which x corresponds. Integrating fllSD and taking 
its expectation, we have 

E[H t ] = H + \k 2 f E \v 2 H s - f3V 8 ] ds. (16) 
2 Jo 1 J 

Let Ut denote E[H t ], the unconditional energy expectation, which can be interpreted as the 
total system energy at time t. Then by differentiating (|I6|) we obtain 

dU t 1 2 



-K 2 ((n + l)(H-U t )-(3E[V t }) , (17) 



dt 

by use of the Laplace-Beltrami equation for H t . This relation can be integrated to yield the 
solution for the time development of the system energy: 

U t = H+(H -H)e-^ 2(n+1)t 

-\k 2 (3 C [ e^ Hn+l){s - t) V{x)p{x,s)dVds. (18) 
2 jo Jr 

In the limit t — > oo the only contribution is given by the trace term H and the integral 
of V(x), and the resulting energy approaches the internal energy of the system in thermal 
equilibrium. In particular, in the high temperature limit (3 — > the contribution from 
V(x) vanishes, and we recover the uniform average of the energy eigenvalues. In the low 
temperature limit (3 —>■ oo, the gradient term in the drift of (0) dominates, and the system 
is forced to fall to the ground state. It is interesting to observe, as a consequence of (|l7|). 
that once thermal equilibrium is reached, we have the identity 

kT 2 C = Var[F] + (n + l)kT(U - H) , (19) 

for the heat capacity C = dU / dT, where (3 = 1/ kT and Var [H] = (AH) 2 is the unconditional 
energy variance. This follows from the conditional variance formula 
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Var[#] = E\yar x [H\]+Vax[E x [H]] , 



(20) 



where E X [H] and Va.r x [H] are, respectively, the conditional expectation and the conditional 
variance of the energy, when the system is in the pure state x. This relation expresses 
the total energy uncertainty by the sum of terms corresponding to quantum and thermal 
uncertainties. 

As an illustration, consider the case of a spin one-half particle in heat bath. The state 
space is CP 1 , which we view as a 2-sphere, the north and the south poles corresponding to 
the upper and lower energy eigenstates, with energies +h and —h, where h is the magnetic 
moment of the particle times the external field strength. The symplectic flow gives rise to 
latitudinal circular orbits on the sphere, while the gradient flow is in the direction along the 
great circles passing through the two eigenstates, pushing the state towards the south pole. 
The equilibrium energy latitude is obtained by balancing the gradient flow and the Brownian 
fluctuations. If we let 9 denote the angular coordinate for the state as measured from the 
north pole, then the squared energy uncertainty, given that the system is in a pure state 
at latitude 9, is Vare[H] = h 2 sin 2 9. The conditional energy expectation is Eg[H] = hcos9. 
Since H vanishes, the only contribution to the energy in flT8|) is from the volume integral, 
which gives U = f3~ x — hcoth(/3h), in agreement with the result of a direct calculation of 
^[^e [//"]]. At infinite temperature the system energy corresponds to that of the equator, 
whereas for finite temperature the system energy corresponds to that of a latitude in the 
southern hemisphere. At zero temperature, the system collapses to the ground state. 



V. ENTROPY PRODUCTION 

Next we consider the dynamics of the total entropy of the system, given by St = 
— J r p(x, t) In p(x, t)dV. Without loss of generality we can write 

n(rf s exp(-/3(H(x)+r](x,t))) 

^ , J f r exp(-p(H(x)+7i(x,t)))dV ' K ' 

for a general density function, where rj(x, t) determines the departure from thermal equilib- 
rium. Substituting this expression into the formula for St, and making use of the Fokker- 
Planck equation (|TTJ) with drift as in (0), we find 



dS t 



at 

On the other hand, for the total energy Ut we have 

dUt 
dt 

Therefore, we obtain 

dSt_ 
dt 



]^k 2 (3 2 j r (V a V V a H + V a V V aV )pdV . (22) 



= \k 2 (3 J r (V a V V a H)pdV . (23) 

rfll 1 r 
P-fr + iff? J r (y a V V aV )pdV . (24) 



This is the general formula for the entropy production associated with the process ( fjj ) 
which shows that 
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f >- ^ 

in all circumstances. In particular, if U t increases as a consequence of heating, then the 
entropy of the system also necessarily increases. 



VI. EVOLUTION OF THE DENSITY MATRIX 

Finally we derive the dynamics of the density matrix Pp(t) associated with the process 
(|H|). Let us revert to homogeneous coordinates for the state space CP n , and write LT^(x) = 
ip a (x)ipi3(x)/ip' r (x)ip 1 (x) for the projection operator onto the pure state x G r represented 
by the Hilbert space vector if) a (x). From equations ([7]) and (|T0|) it follows that for a general 
(possibly nonlinear) observable 4>{x) on F we have 



d_ 

Of 



E[<f>(x t )} = J r {2n ab V a <pV b H + ^K 2 V\ 



~K 2 {3g ab V a <pV b H)rA^i)<iV. (20) 



However, if <p{x) is given, more specifically, by the conditional expectation of an ordinary 
linear quantum mechanical observable 0^, then <p(x) = 0„n^ (x). Thus for a linear observable 
we have E[4>{x t )] = 4>aP^i^)y where 

p%t) = J r U a p (x)p(x,t)dV (27) 

is the time dependent density matrix associated with the state p(x, t). To obtain the equation 
of motion for Pp(t) we need to evaluate the three terms in the integrand on the right of (p6|). 
These are given as follows. For the commutator term we have 

n ab V a <J>V b H = \i{^Hl - Hfyl) , (28) 
for the term involving the Laplace-Beltrami operator we find 

W = - (n + , (29) 

and for the anticommutator term we obtain 



g ah V a <pV b H = -(tfHZ + HfyDUf - (<£H$)I§IQ . (30) 



Inserting these expressions into fl26|), and noting that the result must hold for arbitrary 
we deduce that 



0_ 



p% = i H«p} - p«H} - —k /3 (H«p} + p«H] 



1 



2 •( ( £jr<> , * „i» z^"7 



(SI - (n + l)p a + PH«p«]) . (31) 
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Here the density matrix p%(t) is given by (p7|), the first moment of the projection operator 
Up(x), whereas p^{t) is the second moment of ITg (x), given by 

p£(f) = J r U^x)U](x)p(x,t)dV . (32) 

Equation (^l|) is the dynamical law for the density operator, which takes an arbitrary initial 
state into an equilibrium state, with heat bath temperature T. The first term on the 
right of (0) leads to the Liouville equation of linear quantum dynamics, while the general 
nonequilibrium process has a richer structure. The emergence of the second moment term 



in (31) can also be interpreted by analogy with the renormalisation group equations. 



VII. DISCUSSION 

The model we have described here is surprisingly tractable, and has many attractive 
features. Experimental support could be pursued in two stages. First, we point out that 
it is not difficult, in the case of the canonical T-ensemble, to derive explicit formulae for 
the partition function Z(j3), the state density tt(E), the density matrix p^, and the second 
moment p^]. If these results turn out to give a better account of equilibrium phenomena 
than the conventional approach in some context, then the next step would be to look for 
effects resulting from the higher moment term in the nonequilibrium dynamical law fl3~T|). 
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